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Introduction 


NASA's AES is pioneering new approaches for 
future human missions beyond Earth orbit. 

All spacecraft systems must be minimized with 
respect to mass, power, and volume. 

Here, we show work related to improving system 
efficiency and reliability for water separation 
systems on crewed vehicles 

— the initial development of COMSOL simulations in 
support of the Atmosphere Revitalization Recovery 
and Environmental Monitoring (ARREM) project 



Isothermal Bulk Desiccant 


Simultaneous desorbing and adsorbing flows 

Two paths are thermally linked 

- Adsorption more efficient due to heat loss 

- Desorption more efficient due to heat gain 

Linking can be done in variety of ways 

- e.g., Meshes or grids filled with sorbent 

- Here, simply a thermally conductive housing (Al) 

Flow paths alternate every 'half-cycle' 

Can be driven by different means 

- e.g., vacuum 

- Here, driven by slight pressure differences from 
atmospheric 



Isotherma Bu k Desiccant 


• Geometry choice 
impacts mass 

• Tested a sub-scale 4-cell 
IBD version 

• Used 13X sorbent (large 
and fragile pellets) 

• Both a plastic (unlinked) 
and Al 6061 (linked) 
housing 

• Covered inlet and exit 
with mesh-covered 
baffles 



Isothermal Bulk Desiccant 


Imported IGS CAD file 
directly into COMSOL 

Meshed immediately 

- Working with an un- 
converged mesh 

Red/blue domains are 
paired regions of 
wet/dry flow 

Baffles and manifold not 
modeled here 



The size of the IBD bed in the three 
dimensions are shown in inches 



IBD Comsol Model Description 

Physics Nodes: 

- Free and Porous Media Flow with Forchheimer Drag 

- Dilute Species Transport with Convection, Diffusion, and Reactions 

- Heat Transfer in Porous Media with a Heat Source 

- Heat Transfer in Solids 

- Domain ODE 

- Wall Distance (modified Eikonal equation) 

Transport of a dilute species, H 2 0, in a carrier gas, N 2 , modeled 
free and porous media flow with pressure boundaries. 

The sorption rates and pellet loading were determined from 
solving a distributed ODE. 

Heat transfer in the porous media and the solid housing were 
modeled with the respective physics nodes. 

A wall-distance calculation was required to determine the local 
packing density of the pellets. 

Boundary conditions (inflow/outflow and wet/dry air) were 
switched after every half-cycle of the desiccant bed. 



IBD Comsol Model Inputs 

Porosity model: 


S (D) = £ 0 (1 + (— - V)e~ 4D/dp 

Sr 


’0 

where D=distance to nearest wall, d p =pellet diameter, e 0 =center-line 
porosity 2 3 

Permeability model (BKZ): a :(s) = — 

150 ( 1 -^ 2 ) 

Sorbent mass balance (LDF approximation): dq 


dt 


= k m (q-q) 


where k m = mass transfer coefficient, q*=equilibrium loading 
Reactions: 
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Heat source: Q = -(1-s) — 5H where 8H=heat of adsorption 

dt 
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Forchheimer Drag Coefficient: 



Equilibrium Loading 

• Loading from Toth equations: 
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where a 0 , b 0 , c 0 , and E are the Toth constants for 
H 2 0 on 13X 



Physical Inputs 

Empirically determined: 

- mass transfer coefficient: 0.033 s 1 

- dispersion coefficient: 0.0012 m 2 /s 

- center-line porosity: 0.3 

Physical details of Sorbent and Housing: 

- density 

- heat capacity 

- heat transfer coefficient 

- thermal conduction 

Flow details: flow rates, inlet temperatures, 
inlet and exit pressures, inlet H 2 0 concentration 



Results: Thermal Characterization 


Experiment not 
initially intended 
for this purpose 

Match to data 
requires 'unreal' 
physical 
parameters 

Test article 
boundary 
conditions very 
complex (thermal 
leaks and sinks) 

Nominally 

insulated 



Nominal physical values (k, k, cp, h) result in 
a much too rapid rise (<.2 hrs) to the inlet 
115 °F value. Plastic housing. 


Results: 


Not yet modeled 
with Comsol due to 
time-step crashes 

Test not well 
diagnosed 

Packing not well 
constrained 

Likely suffering from 
dusting and 
clogging 


Breakthrough 
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Results: Time 

Required checks to ensure 
positive quantities 
(concentration, gas density, 
flow rate, pressures, loading) 

13X isotherm is flat at low 
vapor pressures, so large 
initial loading required 

Takes extremely long time to 
run-time converge (i.e. 
months of calendar time!) 

Aluminum housing. 

End of 8 hour experiment. 

Reached symmetry around the J 
ambient temperature. * 

Inlet dew point 8 °C, so peak 
efficiency here is ~93%. 


Dependent 
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Results: Time-Dependent 

• System is 
'unloading' very 
slowly from 
saturated initial 
conditions 


Concentration still 
too high as well 
(only ~60% 
efficiency) 

Using MATLAB 
script to automate 
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Summary 

Just starting to use Comsol for NASA AES 
ARREM problems 

Initial tests with silica gel seem more 
numerically robust 

Very sensitive to initial conditions (e.g., 
loading) and boundary conditions (e.g., 
volumetric flow inlet does not converge well) 

Not yet spatially converged (takes weeks to run 
on a server) 



